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Motivated by the unabating interest in the spin- 1/2 Heisenberg antiferromagnetic model on the 
Kagome lattice, we investigate the energetics of projected Schwinger boson (SB) wave functions in 
the J1-J2 model with antiferromagnetic J2 coupling. Our variational Monte Carlo results show that 
Sachdev's Qi = Q2 SB ansatz has a lower energy than the Dirac spin liquid for J2 > O.O8J1 and 
the q = Jastrow type magnetically ordered state. This work demonstrates that the projected SB 
wave functions can be tested on the same footing as their fermionic counterparts. 

PACS numbers: 



The Heisenberg antiferromagnet on the kagome lattice 
has long been anticipated to realize a spin liquid (SL). Re- 
cently, Herbertsmithite ZnCu3(OH)6Cl2r'"— which con- 
tains kagome layers of spin-1/2 moments, has emerged 
as an experimental candidate with no sign of any order- 
ing down to 50 mK. Interest in spin liquids on the kagome 
lattice has been re-ignited by recent works;^"— where Den- 
sity Matrix Renormalization Group (DMRG) studies find 
a spin disordered ground state with a small gap. A vast 
review of earlier literature is revisited in Ref. 0- Pre- 
liminary DMRG data suggests that the system moves 
deeper into the spin liquid phase upon adding small an- 
tiferromagnetic second-neighbor J2 coupling;^ and this 
supports earlier Exact Diagonalization (ED) study that 
found an increase in the gap for J2 up to 0.1 iS 

In an early study of the nearest-neighbor model us- 
ing large-A^ treatment of Schwinger boson (SB) slave 
particles, Sachde"\»iS found that condensation of spinous 
gives rise to magnetically ordered ground states for spin 
S > 0.26. A subsequent variational study by Sindzin- 
gre et using Resonating Valence Bond (RVB) wave 
functions interpolating between spin liquids and magnet- 
ically ordered states instead found that the former have 
lower energies. However, the spin correlations beyond 
first-neighbors do not agree well with ED results liiii^ 
More recently, the trial energy of the projected Dirac spin 
liquid constructed using fermionic slave particle approach 
was found to lie very close to the ED ground state en- 
ergy in the nearest-neighbor Heisenberg modelr^^ and a 
very recent study extended this to the presence of second- 
neighbor coupling J2i^ However, the observation of an 
energy gap in the DMRG and ED studies suggests that 
the gapless Dirac spin liquid may lose some ground in the 
presence of J2^^ 

To investigate this possibility, we study the energetics 
of a class of projected Schwinger boson wave functions in 
the J1-J2 Heisenberg model with Hamiltonian 



H — Ji ^ Si • Sj + J2 ^ Si • Sj, 



(1) 



where {i,j} and ((«, j)) denotes first and second-neighbor 
pairs. We set Ji = 1 as the unit for energy and consider 





(a) q = SB ansatz 



(b) ^/3 X ^/3 SB ansatz 



FIG. 1: The SB ansatze {A,j} from Ref. [13. The unit cell 
is shaded for each ansatz. All equidistant Aij have identical 
magnitudes, and Aij is positive if an arrow points from site i 
to j. For the q = ansatz, we extend Aij to include second- 
neighbor pairing. The y/3 x ^3 ansatz has poorer energy for 
J2 > 0. 



only antiferromagnetic J2 > 0. The Schwinger-boson 
representation of a spin S is given hy^ 

S-i^6t^,,,fo,,, K^Y.^%^2S, (2) 



where is a bosonic operator, cr arc Pauli matrices and 
K is the number of bosons per site. The Hamiltonian be- 
comes quartic in the bosonic operators, and upon mean- 
field decouplingiSiiSiil leads to the following 



m.t. = 



^E(a^IA+h.c.)+e 



(3) 
(4) 



We treat the "pairing amplitudes" Aij and "chemical 
potential" /i as variational parameters. Equations ([4]) 
are self-consistency relations in the mean field. Fig- 
ure [T] shows Sachdev's ansatze for Aij which have good 
mean field energies in the Ji-only modelJ^ The nearest- 
neighbor Aij are real, and their signs are positive if an 
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FIG. 2: Mean field "phase diagram" for the q = SB 
ansatz.— The phase boundary separates spin liquid from the 
magnetically ordered phases. A spin liquid regime is present 
for physical spin- 1/2 systems with k—1. The new magnetic 
order is complex and is not relevant for this paper. 



arrow points from site i to j. Following Sindzingre et a/,— 
we label these as q = and -s/S x -s/S Schwinger boson 
ansatze, which correspond to Sachdev's Qi ~ Q2 and 
Qi = --Q2 respectivel y (also called [ttHcx, ORhom] 



and [OHex, ORhom] in Refs 
ansatz shown in Fig. 1(a) 



[THi). For the q = SB 

we introduce an additional 
real parameter for the second-neighbor pairing, with the 
pattern of arrows going clockwise in triangular loops for 
both first and second neighbors. Second-neighbor Aij for 
the y/S X -s/S SB ansatz are forbidden by its projective 
symmetry group (PSG). In Ref. [13, Wang et al found 
two more distinct ansatze for symmetric spin liquids, but 
they argued that the Heisenberg model energies of those 
spin liquids are expected to be considerably poorer. Our 
variational calculations confirm that this is indeed true. 

We first present the results of a crude study on the 
accessibility of Schwinger boson spin liquids at the mean 
field level, by computing the critical boson density 



Kc = -1 
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accessible in the mean field without Bose condensation in 
the thermodynamic limit. Here {Aq,} are eigenvalues of 
the matrix —iA [cf. Eq. (|6]) below], /J,inax = — niax{|AQ,|}, 
and N is the number of sites on the lattice. Below Kc, 
the mean field excitation spectrum is gapped and gives 
rise to a stable spin liquid. For k > Kc, the gap closes 
and magnetic ordering results from spinon condensation. 



With only nearest-neighbor Ai 



0.5 and 0.54 for 



the q = and a/S x a/S SB ansatze respectively .^"'^'' In 
analogy to Wang's analysis for the honeycomb lattice^ 
Fig. [2] shows the critical boson density Kc versus A2/A1 
for the q = SB ansatz, where Ai and A^ are the am- 
plitudes of first and second-neighbor Aij . We note that 
Kc > 1 in the parameter range —0.4 < A2/A1 < —0.18, 



i.e., the second- neighbor pairing has opened up a disor- 
dered regime relevant for S = 1/2^ 

We now turn to a variational Monte Carlo (VMC) 
study of the J1-J2 model on the symmetric 36-site clus- 
ter used in previous numerical studieS )^^'^^ which allows 
a direct comparison with ED energies as well as the ener- 
gies of the Dirac SL and magnetically ordered states. We 
construct projected SB wave functions as follows. Writ- 
ing the real anti-symmetric matrix A [see Eq. (|3])] as 



A = iMkM\ 



(6) 



where A is diagonal and M is unitary, we solve the mean 
field Hamiltonian using Bogoliubov's transformation and 
obtain the following trial wave functioni^ 



t ^t 



I^sb) = Vg ey.^ l^ujk h]^h 
Uik = I > — , 



(7) 
(8) 



The Gutzwiller operator Vg enforces the constraint k — 1 
at every site. Although this density may not be accessi- 
ble to a given ansatz and fj, at the mean field level (see 
Fig. the projected SB wave function is a valid vari- 
ational state in the physical Hilbert space. Here, Ujk 
decays exponentially for /i < /imax and the projected 
SB wave function realizes a short-range RVB state when 
written in the valence bond basis. In the limit ^ <^ fJ-max, 
Eq. ^ shows that the pattern of Ujk roughly follows 
that of Ajk- More generally, the PSG of the q = SB 
ansatz enforces the pattern of postulated in Sindzin- 
gre et alr^ where Ujk only connect sites on the different 
"sublattices" A,B,C defined in the sense of the q = 
magnetic order in Fig.ISja). At /i = /imax; Ujk decays in 
a power-law ~ \rj — r^j"'^; one can view /i /^max as a 
finite size realization of magnetic ordersi ^^'^^ 

In our VMC simulations, the amplitude of each sam- 
pled spin configuration is given by the permanent of the 
A^/2 x A'^/2 matrix {ujk}, where j and k run over the 
spin-up and spin-down sites respectively. We make use of 
Ryser-Nijenhuis-Wilf dense permanent algorithm to cal- 
culate the permanent)^ in this way, the procedure does 
not have the sign problem encountered in the valence 
bond basisiii Despite a poor computational cost scaling 
^ 2^/^ with system size, it is manageable for sizes that 
arc already interesting. 

We also consider magnetically ordered states shown 
in Fig. [3] which arise from the condensation of spinous in 
the respective ansatzi^^ For both orderings, their classical 
nearest-neighbor energies are identical, but the second- 
neighbor energy is clearly lower for the q = ordered 
state since it has antiferromagnetic second-neighbor cor- 
relations while the x ^/S state has ferromagnetic cor- 
relations. It is therefore sufficient to consider only the 
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(a) q = MO 




(h) VSx MO 



FIG. 3: Magnetic orderings (MO) which arise from spinon 
condensation in q = and ^/3 x ^/3 SB ansatze. A, B, and C 
are the 120° antiferromagnetic spin orientations. For J2 > 0, 
the X y/3 MO has poorer energy than the q = MO. 
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former. We construct the following trial wave function 

°o) = exp |z ^ 0,5; - K,,S!S] I , (9) 

where (pj = {0, ±27r/3} are the phase angles on the three 
sublattices in Fig. ^a), and Kij are two-body pseudo- 
potentials for the Jastrow factor. Such Jastrow- type 
wave functions are widely used in VMC studies due to 
their simplicity. We allow two variational parameters for 
the first and second-neighbor pseudo-potentials, and two 
more for a power-law decay between further neighbors. 
We also consider a Huse-Elser— type of three-site phase 
factor allowed by the symmetry of the classical state, but 
it apparently does not improve the trial energy. 

To get an idea of the variational energetics landscape, 
we also include the Dirac SL constructed from fermionic 
spinous hopping with flux tt through hexagons and flux 
through elementary triangles J^ii^ This state was ex- 
tended in Ref. [H to include second-neighbor hopping 
such that triangles formed by two nearest-neighbor bonds 
and one second-neighbor bond have flux n. The ampli- 
tude of the second-neighbor hopping provides a single 
variational parameter. For large sizes, we reproduce re- 
sults in Ref. [l5| for J2 > 0; we perform VMC for the 
present 36-site cluster and find negligible size dependence 
on the scale in Fig. 

Figure m shows the variational energies of the Dirac SL, 
the q = SB wave function, and the q = magnetically 
ordered state. For the Ji-only model, the Dirac SL has 
significantly better energy than the SB state. However, 
the latter improves quickly with J2 > 0, and becomes 
lowest for J2 > 0.08 among the wave functions in this 
study. The q = Jastrow-type MO state has higher en- 
ergy for all J2 values shown. 

The q = SB wave function has two variational pa- 
rameters, A2 and /i. For J2 = 0, the lowest energy of 
—0.420 per site occurs at A2 « —0.15. Interestingly, this 
A2 approaches the spin liquid window in Fig. [5] where 
we observe fairly narrow spinon bands. We find antifer- 
romagnetic correlations between second-neighbor sites. 



FIG. 4: Comparison of trial energies per site for Dirac SL, 
q = SB wave function, and q = Jastrow-type magnet- 
ically ordered (MO) state. The SB state has poorer energy 
than Dirac SL for J2/J1 < 0.08, but performs better for larger 
J2 and better than the Jastrow-type MO for all J2. 
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FIG. 5: Optimal A2/A1 versus J2/J1 for the q = projected 
SB wave function, also optimized over jj, for each J2. 

These results are very close to those obtained by Sindz- 
ingre et al in a variational study of the Ji-only modeljii 
wherein they considered RVB ansatze with a few varia- 
tional parameters for nearby-neighbors Ujk and a power- 
law decay ~ \rj — rk\~^ with p = 5 for further neighbors. 
For our projected state, optimizes very close to /^max 
which corresponds to formal p = 3 in the thermodynamic 
limit. Despite the difference in the details of the real- 
izations, both are suggestive of a near-critical state at 
J2 = 0. Thus, Table V in Ref. [Tl| indicates that such 
wave functions have significant q = correlations across 
the full 36-site cluster. 

Figure [5] shows the optimal A2 against the second- 
neighbor coupling J2- We find that A2 increases with 
J2 and is important for improving the trial energy of the 
q = SB wave function. Beyond J2 ^ 0.1, the optimal 
/i starts to decrease away from /Ltmax, e.g., it is 1.02/Xinax 
for J2 = 0.2 and moves further to 1.05/.tmax for J2 = 0.4. 
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Discussion. Our energetics study reveals the q = 
SB wave function as a viable candidate for the J1-J2 
Heisenberg model. This is perhaps not surprising since 
this state is quite competitive in the Ji-only modelii 
and has antiferromagnetic second-neighbor correlations 
which are favorable when J2 > is added. The q = 
SB state can furthermore accommodate the J2 coupling 
by varying A2. In the large J2 limit, the system breaks 
into three independent Kagome networks, each as dif- 
ficult as the original nearest-neighbor Kagome problem. 
The large-^2 state in the large J2 limit is just like the Ai- 
only state in the Ji-only model, so while not the best, is 
again reasonably good in energy. Thus, the A1-A2 ansatz 
provides a way to interpolate between the small J2 and 
large J2 regimes and is an appealing candidate. It wins 
against the Dirac SL and against the best Jastrow wave 
function for the q = MO, but it may also correspond 
to possible MO at intermediate J2- It would be very in- 
teresting to check our results against exact calculations 
in the J1-J2 model on the 36-site cluster to assess the 
accuracy of the projected q = SB state. 

From our Schwinger boson wave function study, we 
cannot address the question whether the ground state 
is spin liquid or has magnetic long range order. It is 
knowniiii^ that the RVB wave functions can realize both 
phases depending on the range of the valence bond am- 
plitudes Ujk- In the projected SB wave function setup, if 
fi is very close to /.tmax, this can be viewed as a finite-size 
realization of the spinon condensation and hence mag- 
netic order. On the other hand, if /.t is a finite distance 
away from /imax, this gives exponentially decaying ujk 
and hence short-range RVB spin liquid. While we find 
that /i optimizes away from /Xmax for J2 > 0.1, these 
small-size results cannot be used to establish the long- 



distance behavior, and the ultimate phase determination 
must come from exact studies on larger systems. Never- 
theless, we hope that our demonstration of the viability 
of the q = SB wave functio n^^i^^i^^'^"^ can be useful for 
further studies of the J1-J2 Kagome antiferromagnet. 

In this work we used a dense permanent routine}^ such 
calculations can in principle be pursued to 48 sites. If we 
also restrict Ujk to only few nearby-neighbors, the VMC 
can be scaled further due to sparseness of the matrix. 
Simulations in the valence bond basis may reach larger 
sizes attention to the sign problem is needed there al- 
though it is less severe than the sign problems in QMC. 
An important aspect of our work is the demonstration 
that projected Schwinger boson wave functions can be 
tested on the same footing as the slave fermion spin liq- 
uids, for smaller but still reasonable system sizes, and 
can be included in the VMC toolbox. Here we high- 
light our use of permanents in a variational study of 
Heisenberg model on the triangular lattice in magnetic 
Heldr^ where we obtained excellent wave functions for 
Mott insulators and supersolids of bosons with frustrated 
hopping. We suggest the honeycomb spin liquid^ and 
Wang's proposal^! as one context for applying the pro- 
jected SB wave functions, as well as other model pro- 
posals in Ref. [l3 for realizing new spin liquids on the 
triangular and kagome lattices 
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